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FLIC Overlap Fermions 
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The action of the overlap-Dirac operator on a vector is typically implemented indirectly through a multi-shift 
conjugate gradient solver. The compute-time required depends upon the condition number, k, of the matrix that 
is used as the overlap kernel. While the Wilson action is typically used as the overlap kernel, the FLIC (Fat 
Link Irrelevant Clover) action has an improved condition number and provides up to a factor of two speedup in 
evaluating the overlap action. We summarize recent progress on the use of FLIC overlap fermions. 



1. Introduction 

Overlap fermions Q are a realisation of chiral 
symmetry on the lattice. Given some reasonable 
Hermitian-Dirac operator H, we can deform H 
into a chiral action through the overlap formal- 
ism, 



Do = l{l + 75e{H)), e{H) = 



H 



(1) 



Unfortunately, the matrix sign function e{H) is 
difficult to evaluate and is typically approximated 
by a sum over poles |^ which can be evaluated us- 
ing a multi-shift conjugate gradient (CO) solver 
|3j. This is an iterative approximation where 
the number of iterations for a given accuracy in- 
creases with the condition number of the kernel, 
i^{H) = I Aiiiax/Aminl- 
Usually the Hermitian Wilson-Dirac operator is 
used as the overlap kernel. Its low-lying spectrum 
is characterised by a handful of isolated eigen- 
modes which can be very small, increasing the 
condition number, k, unacceptably. These eigen- 
modes can be projected out of the basic operator, 
reducing its condition number to a numerically 
acceptable level, and then dealt with explicitly 
m. Unfortunately, as the spectrum rapidly be- 
comes dense, projecting out low-lying modes can 
only help one so far. 

An alternative is to use a kernel with an im- 
proved spectrum, that is, where the region of 
dense modes is shifted away from zero. The FLIC 
(Fat Link Irrelevant Clover) ||J^ action possesses 



this property. 



2. Spectral Flow and Condition Number 

Spectral flow diagrams give us a good compari- 
son of the two different kernels, FLIC and Wilson, 
as they allow us to see directly the difference in 
the qualitative structure of the low-lying spectra 
of the two actions |^ . Figure |l| displays the flow 
of the lowest 15 eigenvalues as a function of m for 
an ensemble of 10 mean-field improved Symanzik 
configurations at /3 = 4.60 and size 12^ x 24, with 
a = 0.122(2). As we are interested in the mag- 
nitude of the low-lying values rather than their 
sign, we illustrate |A| vs m. 

We see the Wilson spectrum is very poor, with 
a high density of very small eigenmodes and no 
gap away from zero. By contrast, the FLIC ac- 
tion (which has a clover term, and irrelevant oper- 
ators constructed from four-sweep APE-smeared 
links) possesses a spectrum which is clearly su- 
perior to that of the Wilson. Not only has the 
density of very small modes been significantly re- 
duced, the region where the spectrum becomes 
dense has been shifted much further away from 
zero. 

This improvement in the spectrum results in 
the FLIC action being much better conditioned 
than the Wilson action. In Figure || we exam- 
ine the condition number, /t, of the two actions 
as a function of m. The condition number is cal- 
culated after having projected out the lowest 15 
eigenmodes on the 2 lattices that are shown. The 
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Figure 1. Spectral flow of the Wilson action (left) 
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Figure 2. Comparison of the condition number 
improved glue (left plot) at both a = 0.165 and a - 
(right plot) at a = 0.165 with 10 modes projected 
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and the FLIC action (right) at /3 = 4.60. 
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the FLIC and Wilson fermion actions. Symanzik 
0.122 with 15 modes projected out, and DBW2 glue 
ut are illustrated. 



finer lattice is the same as the one used in the 
spectral flow plots, and the coarser lattice uses 
the same gluonic action, but is an 8"^ x 16 lattice 
at /? = 4.38, corresponding to a lattice spacing of 
a=0. 165(2). The points are the mean condition 
numbers across the ensembles, and the error bars 
indicate the minimum and maximum condition 
numbers, giving an idea of the variation in k. 

3. Gluonic Action 

There have been suggestions to accelerate the 
computation of the sign function by using non- 
pertubatively improved gauge actions ||^. Our 



results are based on Symanzik improved gauge 
configurations and further improvements arising 
from the use of the FLIC action are in addition 
to that of using improved gluon-field configura- 
tions. This is verified by performing a similar 
analysis on gluonic configurations using a Monte- 
Carlo Renormalisation Group improved action. 
Some preliminary results of this investigation us- 
ing DBW2 glue § are displayed in Figure ^, with 
a full report given elsewhere We note that at 
a fixed lattice spacing, DBW2 glue improves the 
condition number for both actions, although the 
effect is much more pronounced for the Wilson 
action. 
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4. Compute Time 

Saving iterations (by reducing the condition 
number) does not necessarily reduce the most im- 
portant quantity, compute time. Shifting from a 
standard Wilson action to a partially smeared ac- 
tion means that we now have two sets of gauge 
fields, the standard and smeared links. Ad- 
ditionally, the standard spin-projection trick is 
no longer applicable, possibly providing an ad- 
ditional factor of two in compute time needed. 
However, it can be shown that the spin projec- 
tion trick can be generalised to include partially 
smeared actions as well |10 . This results in pay- 
ing at most a single factor of two compute-time 
for a FLIC-fermion matrix-vector multiplication. 
As there is significant additional expense in the 
evaluation of the overlap sign function, we get to 
keep the majority of the speedup gained by re- 
ducing the number of iterations This results 
in FLIC-Overlap fermions being approximately 
twice as fast as the standard Wilson-kernel for- 
mulation. 

5. Physical Results 

Regardless of the kernel used, all overlap 
fermions are free from 0{a) errors. However, 
different kernels may in general produce actions 
which differ at O(a^). As a first investigation into 
this matter, we have calculated the quark prop- 
agator in momentum space using FLIC-Overlap 
fermions, essentially performing the same calcu- 
lation that has been done earlier with the stan- 
dard overlap action The chirally extrapo- 
lated mass function is shown in Figure ^, with a 
full report given elsewhere |12 . 



6. Conclusion 

In the overlap formalism one is free to choose 
the argument of the sign function, the overlap 
kernel, so long as one uses a reasonable Dirac op- 
erator. The standard kernel choice is the Wil- 
son action. By choosing the FLIC action, one 
can obtain a significant reduction in the compute 
time needed for overlap fermions. We have re- 
viewed the work that has been done to date us- 
ing FLIC Overlap fermions, including some pre- 
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Figure 3. The chirally extrapolated mass function 
as a function of p, at P — 4.60, L — 12^ x 24. 

liminary results into the physical structure of the 
quark propagator. 
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